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ABSTRACT 


Waves in the surf zone are a highly nonlinear process 
which is evident by the appearance of secondary waves. The 
secondary waves appear as strong peaks in the period PDFs 
corresponding to the first harmonic of the peak of the wave 
spectrum. The strong first harmonic period peak is also 
Goel cred i nmene Milach Ivscorrelared height and velocity PDis. 

Preaking waves may be viewed as a gradation between two 
extremes. The period, height and velocity PDFs for spilling 
breakers generally were unimodal, whereas for plunging 
breakers strong bimodality was found. 

The joint probability density functions for periods and 
heights of the breaking waves show high correlation (0.60- 
0.80) which says that greater wave periods are associated 
with larger breaker heights. The joint PDFs of period and 
particle velocity, and velocity and height, suggest that the 
maximum onshore particle velocities are correlated with both 
the wave periods and wave heights. 

Due to the high probability of the secondary waves, the 
mean wave period for breakers is a poor descriptor of the 


average period of the offshore incident waves. 
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I. INTRODUCTION 


Waves in deep water can generally be represented to a 
first approximation as a Gaussian process having zero skew- 
ness and kurtosis equal to three. The random sea surface 
is viewed as the linear superposition of statistically inde- 
pendent spectral components which is adequately described 
by the energy density spectrum. The restriction of statis- 
tical independence implies the wave components act as free 
waves and there is no coupling or interaction between 
various spectral components. 

As the waves begin to shoal, the crests become more 
peaked, and the troughs become more elongated; this results 
in a positive value of the third statistical moment, i.e., 
the probability density distribution of sea surface eleva- 
tions becomes skewed to the right. The skewness is a result 
of phase coupling of the wave components, which is indica- 
tive of the increasing importance of the non-linearities 
or the system. This initially weakly non-linear system 
found in shallow water may be adequately described by a 
Gram-Charlier probability distribution (Longuet-Higgins, 

OG 3). 

Near breaking waves become asymmetrical about the crest 
with a continually steepening of the wave face. Miller 
(1976) has shown that a jet is then formed at the wave 


crest at breaking, which extends forward and then down 





until it closes with the forward face of the breaking wave 
form. The asymmetry about the wave crest yields a change 
in the fourth moment, or kurtosis, which is also indicative 
of the increasing non-linearities. Thornton et al. (1976) 
have noted that secondary wave formation becomes apparent 
on the trailing face of the wave form, due to transfer of 
energy from the primary frequency to its harmonics which 
also points to a strongly nonlinear system. 

This thesis uses data collected as described by Galvin 
(1975) and Olsen (1977) plus additional data collected by 
the author. The experiments were designed to measure the 
sea surface elevation and the onshore particle velocities 
by means of capacitance wave gauges and electromagnetic 
flowmeters, respectively. The experiments were conducted 
at a variety of beaches in order to include wave types of 
spilling, plunging, and collapsing forms. The characteris- 
tics of the breaking wave forms are a function of the off- 
shore wave spectrum and the nearshore bottom slope and 
topography, where the increasing beach slope increases the 
rapidity of the breaking of the waves. 

Empirical distributions were obtained for the periods, 
wave heights, and maximum onshore particle velocities, for 
waves which had been defined by the zero-up-cross method. 
Each of these distributions was compared with the Rayleigh 
probability density function (PDF). Joint PDFs of periods 
and wave heights, maximum velocities and wave heights, and 


perlods and maximum velocities were calculated. 





I~. HISTORICAL BACKGROUND 


Ao OUAIATIVE DESCRIPTION 

Thornton et al. (1976) have described a number of ob- 
served universal similarities for the various types of 
breakers occurring on different beaches. Normally there is 
a quick drawdown of water just before the breaker arrives, 
followed by a steep, nearly vertical leading edge. Miller 
(1976) suggests that subsequent to the development of the 
leading vertical edge, a jet forms at the crest. This jet 
progressively extends itself forward and down wntil it wnites 
with the forward slope of the breaking wave, resulting in 
the formation of a vortex with a circulation around a 
cavity. When the tip of the "jet'' touches the forward slope 
of the breaking wave, the resulting force forms a second 
weaker vortex in front of the first one. In this manner, a 
series of vortices yields a series of large, single pulse 
eddies, which contribute to the rapid decay of the wave 
during breaking. It has been found that the magnitudes of 
breaker vortices are generally correlated with a gradation 
from plunging to spilling breakers. 

Thornton et al. (1976) have observed a sloping profile 
toward the trailing edge of the wave. Secondary waves are 
frequently noted on the trailing edge both on the analog 
strip charts, Figure 1, and by visual observations in the 


field. The secondary waves have been observed to grow during 
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the shoaling process, especially just prior to breaking 
when the primary waves are steepening significantly. The 
secondary waves are a result of the transfer of energy from 
the primary wave frequency to harmonics, and help maintain 
the potential energy across the surf Zone. 

Gallagher (1971) has photographed secondary wave forma- 
thon occurring across reefs. He suggests that the secondary 
waves are brought about by a response to vertical accelera- 
tion near the front of a broken or breaking swell. When the 
amplitude and surface curvature of a shallow water wave ex- 
ceed certain values, the non-linearities of the vertical 
accelerations bring soliton wave forms into being. 

Figure | shows that the maximum particle velocity occurs 
near the time of maximum surface elevation. It can also be 
noted in Figure 1 that there are secondary velocity maxima 
which correlate with the secondary surface elevation maxima. 
The secondary velocity peaks also are a result of the non- 


linearity associated with the breaking wave. 


Bee GIIERGY DENSITY SPECTRA OF BREAKING WAVES 

In order to better understand the PDFs as determined in 
this study, the energy density spectra measured by Galvin 
(1975) and Olsen (1977) were reviewed. The energy density 
spectrum 1s proportional to potential wave energy within a 
frequency band. The energy density spectrum, cca, indicates 
how the variance, a, is distributed over frequencies, f, 


where 


eZ 





co 


Z 
Oo, -| G,(f) df a) 


O 


(Kinsman, 1965). 

The energy density spectra, as estimated from the sea 
surface elevation records, were used to find the frequency 
of peak density and the occurrence of harmonics. The power 
spectra that were measured by Galvin showed a peak density 
@ver tne Erequency range of 0.061 to 0.092 Hz (16.4 to 10.9 
sec). First, second and third harmonics can be observed 
as shown in the spectra of surface elevation and horizontal 
velocities for 4 March 1975. Additionally, subharmonics 
were observed at 0.011 and 0.922 Hz in the velocity spectra. 
Galvin (1975) has eee ees ted that the observed low frequency 
sub-harmonic peaks may be due to edge waves in the surf zone. 

Thornton et al. (1976) note that the appearance of strong 
harmonics could be real due to secondary waves or artificial 
due to Fourier analysis of a peaked wave form. It is noted 
that secondary maxima of n are apparent in the strip chart 
record shown in Figure l. It is felt that the secondary 
spectral peaks are real since secondary waves appear on the 


Sieaip Charts. 


oe OSA ein Y SBENSERY FUNCTIONS OF PERIODS AND RELGHTS 
studies of the PDFs of period, height, and wave leneths 


have been conducted over the last 20 years. The majority 
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of these studies have been mainly concerned with theoretical 
or deep water aspects of the problem. 

The theoretical distribution of wave amplitudes for a 
narrow-banded Gaussian surface elevation in deep water, was 
shown by Longuet-Higgins (1952) to be a Rayleigh distribu- 
tion. Tayfun (1977), in studying the transformation of deep 
water waves to shallow water waves, showed that the Rayleigh 
distribution for wave amplitude was generally applicable to 
all bandwidths. Using the assumption that the wave height 
is twice the wave amplitude, the wave height PDF is then re- 


presented by: 





4 
p(H) dH = = exp Es dH (a) 


where o is the variance of the surface elevation. Using 
pressure records in the Gulf of Mexico, Longuet~-Higgins 
(1975) observed that the Rayleigh distribution fit the ob- 
served distribution reasonably well. He found that there 
is a slight excess of waves with heights near the middle 


On 


the range and a deficit at the two extremes. Since much 
of the high frequency portions of the wave records are fil- 
tered out by the pressure transducer, Longuet-Higgins (1975) 
Suggests that the narrow-band approximation may not be as 
applicable for the wfiltered records. In shallow water 
with much steeper waves, the distribution can again be 


expected to be less applicable due to the non-linearities 


which become increasingly important. 


ies 





Since the statistics of H and T for a broad-~banded 
spectrum were still somewhat unclear, Goda (1970) simulated 
wave profiles using an infinite series where the amplitudes 
were specified by various theoretical spectra and the phase 
was random with a computer. He then examined the simulated 
wave records for surface elevation,crest-to-trough wave 
heights and zero-up-cross wave heights. He found that for 
the zero-up-crossing determination of wave heights, the Ray- 
leigh distribution is a good approximation irrespective of 
the spectral width. 

Goda (1970) suggests that the deviation of the wave 
height distribution from the Rayleigh type in shallow water 
should be considered separately from the change in the func- 
tional shape of the wave spectrum. It must be a highly non- 
linear phenomenon, which requires full analysis with the 
theory of non-linear wave spectra. 

Bretschneider (1959) investigated the PDFs of wave 
heights, H, wave lengths, L, and wave periods, T. He con- 
cluded from his empirical study that in deep water wave 
systems, the PDFs of wave heights and lengths can be approxi- 
mated by a Rayleigh distribution. Using the relationship, 
wave length proportional to period squared for deep water 
linear waves, he concluded that the PDF of periods squared 
is a Rayleigh PDF given by 

2 


4 
p(T) = 2.74 exp |-0.675 A (3) 


where T is the mean period. 
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Using simulated wave data, Goda (1970) found that the 
crest-to-crest periods squared were generally not well ap- 
proximated by the Rayleigh PDF. The mean and standard devia- 
tions of the periods by this method were found to be a 
function of the spectral width. He found that the periods 
determined by the zero-up-cross method are only in fair 
agreement with Bretschneider (1959). If the waves are of 
the type which have sharply peaked spectra, double peaked 
spectra, or a flat spectrum, Bretschneider's (1959) distri- 
bution was totally in disagreement with that of the simu- 
lated waves. 

Contrary to Bretschneider (1959) and Goda (1970), 
Longuet-Higgins (1975) noted that in general the deep water 
wave periods were bell shaped (but not Gaussian). He defines 
the period PDF when calculated by the zero-up-cross method 


as, 


i 


C0) eye 
2(1+n~) 


Sela) are (4) 


No] lo 


where <t> is the mean period and vy is proportional to 


the spectral width where 
= 4 
v= Wy/u) *(<t>/2T) (5) 


and Ho and Hy are the zero and second moments of the 
energy spectrum. 

Koele and de Bruyn (1964), Goda (1967), and Siefert 
(1970) have observed that the distribution of wave heights 


in shallow water does not correlate well with the Rayleigh 


Ih 





distribution. Kuo and Kuo (1974) have suggested that this 
is due to: 1) the non-linear effects of wave interactions 
yielding more larger waves, 2) the ertfect of bottom friction 
yielding a reduction in the low frequency components, and 

3) the effects of wave breaking, which would trumcate the 
distribution and transfer some of the kinetic energy to the 
high frequency components. Therefore the PDFs of wave 
heights with a certain intensity, may be considered to ap- 
proximate a conditional Rayleigh distribution truncated by 
the local breaking height. | 

A zero-up-cross analysis was conducted on waves measured 
on a reef by Black (1978), who then compared the observed 
distribution with the Rayleigh, truncated Rayleigh, and 
the Weibull distributions. He observed that the truncated 
Rayleigh distribution fit well for those waves in the breaker 
zone when the heights were depth limited, The triumcated | 
Rayleigh distribution did not fit well for those waves 
measured either inshore or offshore of the breaker zone. 

He has found that a Weibull distribution more closely approx- 
imated these measured distributions. 

It would be expected that the Weibull distribution would 
yield the better fit. The Rayleigh distribution is a func- 
tion of the variance of the data, whereas the Weibull dis- 
tribution is a function of higher moments about the mean. 
Neither theoretical distribution adequately describes ob- 
served distributions which are multi-modal. 

Black's (1978) observations of wave periods showed a 


distribution which was skewed to the right. This does not 
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fit well with the symmetrical Longuet-Higgins (1975) distri- 
bULlOn, Le was found that the period distribution more 
closely followed the Rayleigh distribution than the period 
Square as proposed by Bretschneider (1959). It was found 
that the Weibull distribution most closely fit the empiri- 
cal PDF, but there is considerable variation in the peaked- 
ness parameter. His attempts in relating the period statis- 
tics to the position on the reef were unsuccessful. He 
feels that this is due to the non-linearities of the waves 
in shallow water. 

In this thesis, the empirical data collected for breaking 


waves were compared with the Rayleigh distribution. 


Paro Roba bicinyY  DENSLTLY FUNCTIONS 

A study of the joint PDF for wave heights, n, and lengths, 
A, which can be viewed as heights and periods) was accom- 
plished by Bretschneider (1959). The general form of the 
joint PDF for two related variables (n,A) is the product of 


ene marzinal PDF (p(n)) and the conditional PDF (p())). 
a 


p(n,aA) = p(n) *p, OQ) (6) 
then 
p(n) = frona = p(n) fe,oa (7) 
© O 


Assuming that the marginal PDFs are Rayleigh, Bret- 
fenne nen eho 7) pexaminaed the Extreme cases of 0 and +1 cor- 


relation. For the case of zero correlation, the marginal 


1s, 





Rayleigh PDFs are given by 


ea) 7 n exp (-1™m*/4) ; n= H/H (8) 
p(A) = BA exp(-1A°/4) 5 A = L/L (9) 


For zero correlation, i.e. total independence (Figure 3), 


the joint PDF is, 
oi e= pica “Pp CAD C0» 


so that the joint PDF of heights and lengths is given as, 
ne EE Z 
p(n,A) = gn exp(-mn /4)° A exp(-TA/4) Cir) 


In terms of heights and periods, equation 15 is, 


3 


p(n,t) = 1.35 Texp(-0n-/4)+* 1” exp(-0.675 t*) (12) 


where t = T/T. 

For the case of total dependence, the correlation coef- 
ficient is equal to one. For a correlation coefficient of 
one, all data points on a plot of joint Rayleigh PDFs fall 
on a 45 degree straight line passing through the origin, 
n=A=0. Assuming that a joint Rayleigh distribution is appli- 


cable, the following relation applies: 
Des Z TT Z 
p(n,A) = 5 n exp(-™n"/4) + y A exp(-7A"/4) Gis) 


which is plotted in Figure 3. 
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Collins (1970) studied the PDFs of breaking waves de- 
rived from deep water waves having a joint Rayleigh distri- 
bution of heights, lengths and angles. He assumed that 
H = yd , where H is the breaker height, d is the breaker 
depth, and y is the breaking criterion assumed. The limit- 
ing factor of waves breaking in the surf zone yields a dis- 
tribution better approximated by a truncated Rayleigh dis- 
tribution. As expected, he found that the longest deep water 
waves were associated with the largest breaking waves. 

Collins (1970) concluded from the study of the joint 
PDFs of wave heights and lengths, that the effect of deep 
water angle of approach is relatively minor for breaking 
waves. The breaking wave steepness was found to be very 
sensitive to the choice of breaking criteria. 

A theoretical expression for the joint distribution of 
wave periods and amplitudes was derived by Longuet-Higgins 
(1975) for a narrow frequency spectrum 


oa SER 
ote iy) = eo sie =e (Lone) 2) (14) 
Zi 





where §€ and n are the normalized amplitudes and periods 


as defined by, 


S= ayy. Gy 
and 

n = (t-<t>)/v<t> CES) 
where v is defined in equation 5. The wave amplitude 





ce 2s is the zero™ 


. moment of the energy spectrum, 


and <t> is the mean wave period. Longuet-Higgins (1975) 
compared this theoretical distribution of ocean waves to 
that of the deep water data obtained by Bretschneider (1959) 
and found good agreement. 

Chakrabarti and Cooley (1977) compared the theoretical 
joint distribution of Longuet-Higgins (1975) with wave data 
recorded during a 1961 storm in the North Atlantic. The 
spectrum of these waves did not fall in the category of 
narrow-band. Through the comparison, it was found that 
there was a considerable agreement at the higher wave heights, 
but a definite trend away from the theoretical curve at the 


lower wave heights. 
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Pil. Exe IMENT 


eee ean NTA SLtES 

Three different Califormia experimental sites were used 
to obtain data on the three types of breaking waves--plung- 
ing, spilling, and collapsing. During the period 4 March to 
10 March 1975, plunging and spilling waves were measured by 
Galvin (1975) on Del Monte Beach, Monterey, Califomia. 

This locale was chosen because the waves are highly refracted 
and directionally filtered, yielding narrowbanded swell off- 
shore which impinges almost perpendicularly on the beach. 

On 29 May 1975 Galvin (1975) conducted a second experi- 
ment on the Carmel River Beach, which is located approximate- 
ly five miles south of Monterey, California. This area is 
in an embayment, where the waves are again highly refracted 
and directionally filtered. The beach at this site is very 
steep resulting in plunging and often collapsing and surging 
breakers. Due to the steepness of the beach and the rapidity 


C= 


breaking of the shoaling waves, reflected waves are found 
to be present. 

Spilling breakers on a gently sloping beach were measured 
MM ascend 7/ eat lornrey -ines Beach, a site just north of 
Paedollas -Califormia (Olsen, 1977). 
'On 83 and 9 March 1978 data were collected by this author 


on Del Monte Beach, Monterey, California. On the first day 


the experiment, very narrowbanded swell was observed 
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breaking on the shore, mainly of the plunging type. A storm 
arrived late on 8 March and continued through 9 March. Winds 
of approximately 15 knots yielded a very confused sea. 
Breakers were of both the spilling and plunging type and 
were observed to break at various small angles relative to 
the beach. Sea surface elevations and velocities were col- 
lected for approximately three hours on each of these days 


on both the flowing and ebbing sides of high tide. 


Bee Loo RUliEN aS 

Capacitance type gauges were used to measure the change 
in the sea surface elevation as explained by Thornton et al. 
(1976). The wave gauges were constructed of 3/8 inch poly- 
propylene covered stainless steel rod. Calibration was 
accomplished in the laboratory statically by lowering the 
wave gauge into a deep water tank. Depth of immersion was 
then plotted against volts. 

The velocity of the water particles was measured with 
Marsh-McBirmmey Model 721 and /22 Electromagnetic Current 
meters. These flowmeters work on Faraday's principle of 
electromagnetic induction and are capable of recording the 
particle velocities in orthogonal directions in a plane. 
The meters were dynamically calibrated in the laboratory. 
Calibration was accomplished by oscillating the current 
meters in a water tank with a prescribed motion simulating 
the flow conditions which would be expected to be found in 
the surf zone. The flowmeter voltage output was found to 


be linear in the velocity ranges found in the breakers. 
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All instruments were placed on towers in the surf zone 
at low tide at the expected point of breaking of the waves 
for the subsequent high tide. Each tower held a capacitance 
wave gauge and an electromagnetic flowmeter. The flowmeter 
was placed below the wave gauge, approximately 0.5 m from 
the bottom, to insure constant submergence. Data were col- 
lected within an hour or two on both the ebb and flow sides 
of high tide. Analog recordings were made on both a strip 


chart and magnetic tape. 


Ge DATA ANALYSIS 

All signals were digitized from their analog form by 
means of a Vidar Corporation 32-channel digital data acqui- 
Sition system. The digitized data were then processed on 
an IBM 360/67 computer. Record lengths of approximately 
30 minutes from each data set were analyzed. Specific de- 
tails of the analysis can be found in Appendix A. 

The data were edited by excluding any value greater 
than six standard deviations from the mean and replacing 
that value with the previous value. The values were then 
linearly detrended to exclude the effects of the rising or 
falling tides. The mean and variance of sea surface eleva- 
tions were calculated for each record as well as the mean 
and variance of the heights, periods and maximum velocities 
all of which are defined on page 27. 

It was desired to only examine the sea-swell band of 
Erequencies lying between 0.03-1.0 Hz (30 sec-l sec). Hence, 
the data were passed through both a high and low pass digital 


ie Leer. 
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The high pass filter had a cut-off frequency of 0.03 uz 
(30 sec). The high pass filter used a Fast Fourier Trans- 
form (FFT) algorithm to obtain the Fourier spectrum of the 
entire 30 minute record. The Fourier coefficients correspond- 
ing to 0-0.03 Hz were then used to synthesize a low frequency 
time series which was subtracted from the wave record. 

The low pass filter was a 25 weight inverse transform 
Filter (Davidson, 1970), as shown in Figure 4. A cut-off 
frequency of 1.2 Hz was used with a terminal frequency of 
gener 

Wave heights were determined from the surface elevation 
record by means of the zero-up-cross method. The zero-up- 
cross technique defines a wave period as the interval between 
adjacent upcrosses. The maximum and minimum of the surface 
elevation within the interval defines the crest and trough 
of a wave. Maximum onshore particle velocities were calcu- 
lated from the flowmeter records. Maximum onshore velocity 
was defined as the maximum velocity amplitude occurring be- 
tween succesSive zero upcrossings as defined by the sea 
surface elevation record. Joint PDFs of wave height and 
period, wave height and maximum velocity, and maximum veloc- 
lity and period were calculated from the analyses. 

Heights and maximum velocities were normalized by divid- 
ing by their respective standard deviations. The periods 
were normalized by dividing by their mean period. The nota- 
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where heights (H'), maximum velocities (V'), and periods 
(T') are the values as determined by the zero-up-cross 
method. T is the mean period, and o, enic ee ceeeaTemetic 
standard deviations of the surface elevations and particle 
velocities respectively. Histograms were calculated for 
the normalized heights, periods, and maximum velocities. 
The histograms were transformed into PDFs by dividing the 
number of value over an interval by the total number of 
values in the array and by dividing by the interval width. 
The mathematical expression for the exact PDF is 

p(x) = lim eT = cia ||; (18) 
where Nee is the number of times the variable occupies an 
interval, N is the total number in the record observed, 
Ax is the interval width and p(x) is the PDF. A restric- 


Eton on the probability density is that 


Joint probability density functions of periods versus 
heights, velocities versus heights and periods versus 
velocities were calculated in a similar fashion. A 15 by 15 
grid was used to determine the frequency of values which 
Simultaneously fill each joint interval. Mathematically 
the probability that the arrays x and y simultaneously 


assume values within the range (x, xtAx), (y, ytAy) may be 


A, 





obtained by taking the ratio of the total number of times of 
eaemeoeccurrence Within a bin, Mey? and the total number of 
observations, N . The exact probability density is approached 
as the number of observations approaches infinity the interval 


widths approach zero. The joint PDF can be defined as: 


lin iL lim Ney - 
2 isi 2 = @ ieee (19) 
yr @ 


Witen the restriction 


-+-co = 00 
fl / px.) dx dy = 1 
O O 


The energy density spectrum (power spectrum) was obtained 
by means of a Fast Fourier Transform (FFT) algorithm. The 
energy density spectrum tells how the variance is distributed 
with respect to frequency. 


The spectral width parameter was calculated using the 


Mpm,-m, 
a a (20) 
04 


where My, My; and m, are the spectral moments (Cartwright 


equation: 


and Longuet-Higgins, 1956). For e=0 the waves have only one 
Frequency and the energy spectrum would be represented by a 

spike. As e+l the sea surface can be described as a broad- 

band process for which the distribution would be Gaussian 


and the spectrum would approach a horizontal line. Swell 
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can therefore be described as a narrow-band process and seas 
as broad-band. Due to the appearance of secondary waves, a 
narrow-banded deepwater spectrum of waves broadens in the 
surf region. 

iemorue laltOnmeoctrrelents On ite  j/Oint values were 
calculated to obtain a measure of the amoumt of linear rela- 
tionship of random pairs of values of periods and heights, 


velocities and heights, and periods and velocities, 


a x 


where x and y are the sample means, x, and y, are 


the sample standard deviations. 


oul 
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Miller (1976) concluded that breaker shapes form a con- 
tinuous gradation of shapes from spilling to plunging. Each 
breaker type undergoes the same processes, Dwe "dil brews Jom ly, 
in the magnitude of the vortices generated. The probability 
density functions appear to show a similar gradation. Two 
simplified cases, for spilling and plunging, were chosen to 
represent extreme examples and are used for all discussions. 
The 4 March 1975 data most simply represent the plunging 
case and the 16 March 1977 data represent the spilling case. 
The results for all data are given in Appendix C. 

In both examples for spilling and plunging breakers, the 
offshore waves appeared as narrow band swell. A spectral 
analysis of the surface elevation records in the surf zone 
yielded spectral widcth parameters ranging from 0.94-0.97 
(see Appendix A), suggesting a broad band spectrum in the 
surf zone. The broad spectral width can be attributed to 
the non-linear wave-wave interactions yielding harmonics in 
the breaking zone. 

In order to determine the origin of the large spectral 
width parameter, ¢« was calculated in the ranges of 0-0.25 
Hz, 0-0.50 Hz, O-1.0 Hz and O-nyquist frequency. The analy- 
sis revealed that large values of the spectral width para- 
meter occurred in all cases. In the case of 0-0.25 Hz which 


contains the first three harmonics, ¢« was greater than 0.7, 
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which was the case with the smallest «. It is concluded 
that the change from a narrow deep water spectrum to a 
broad-band shoaling spectrum is primarily due to generation 


of secondary waves at harmonic frequencies. 


A. PERIOD PROBABILITY DENSITY FUNCTIONS 
Bretschneider (1959) made the assumption that the wave 
lengths were Rayleigh distributed; using this assumption 


and the linear wave theory relationship 


- £ 7% 
L a T . (22) 
the wave periods squared would be Rayleigh distributed in 
deep water. Bretschneider's analysis can be extended to 

the shallow water case using the shallow water approxima- 


tion for wave length 
ea yen 1 Ze) 


Assuming the wave lengths are Rayleigh distributed, the 
periods in shallow water would be expected to be Rayleigh 
distributed. 

The period PDFs for plunging and spilling examples are 
shown in Figures 5 and 6 respectively. The Rayleigh PDF has 
been superimposed on the period PDFs for comparison. The 
formulas used for the Rayleigh PDF for various parameter 
transformations are given in Appendix A. In both the spill- 
ing and plunging wave examples, it is noted that the period 


PDFs only loosely follow the Rayleigh PDF. The Rayleigh PDF 
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overpredicts the wave periods greater than the Rayleigh mode 
value, and underpredicts the observed periods less than the 
peak Rayleigh value. The truncation of the long period 
waves is probably due to the larger period waves breaking 
farther offshore and then reforming to be recorded by the 
Stationary wave gauge as higher frequency waves. 

The spilling example for 16 March 19/77, shown in Figure 5, 
has the mode lying close to the maximum of the superimposed 
Rayleigh PDF. The peak of the breaking wave distribution 
is located at the first harmonic of the offshore wave period, 
which is attributed to the secondary waves. A peak corre- 
sponding to the peak offshore wave period is absent. 

The plunging waves, as represented by the 4 March 1975 


record in Figure 6, are multimodal. Since the Gaussian, Gram- 


Charlier, Rayleigh and Weibull PDFs are unimodal, none of the 


theoretical PDFs will adequately describe multi-modal PDFs 
HOteenis Dakticular plunging sunk comdition. 

Two strong peaks are found to occur in the plunging 
breaker period PDF. The longer period is found to correlate 
well with peak period of the wave spectrum. The shorter wave 
period appears at approximately the first harmonic. Since 
the center period of the bin is used to represent the inter- 
val, the magnitude of the peaks and the proximity to har- 
Monics as represented by the labels is a function of the 
interval width and the interval cutoffs. A 30 bin PDF was 
calculated for better resolution to determine if the peaks 


were harmonics and Figure 7 shows that they are. 
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The waves periods as determined by a zero-up-cross of 
the surface elevation record for spilling conditions yields 
a poor indication of the peak of the energy spectrum which 
corresponds to the mean offshore wave period. In the spill- 
ing breaker case, the mean period corresponded to the first 
harmonic of peak in the energy spectrum. For the plunging 
breakers, the mean period corresponded to a valley in the 
period PDF between the primary and first harmonic periods. 
The appearance of strong peaks at one half the offshore wave 
period is indicative of the importance of non-linearities 


in the form of secondary waves in the surf zone. 


Peek Cia RObABIiGiY DENSELY RUNCTIONS 

Wave height PDFs for the plumging and spilling wave ex- 
amples are shown plotted against the theoretical Rayleigh 
PDF in Figures 5 and 6. The height PDF for spilling waves 
1s not as easily described as the period PDF. The height 
PDF is loosely described by the Rayleigh PDF. The major 
height peak is centered at 2.2. The Rayleigh PDF overesti- 
mates the density of the small wave observations and tends 
to underestimate the density of the large waves. The heights 
are truncated at approximately six standard deviations. Since 
the wave gauge is stationary and therefore the depth of the 
water is essentially Pence dea th? SeERuNecaelom Occurs due fo 
the breaking of the large waves seaward of the wave gauge. 
The calculated mean height value on 16 March 1977 lies at 
2.7, which actually corresponds to a dip in the PDF relative 


to the peak intervals on either side. 
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The mean height observed in the 4 March 1975 plunging 
breaker PDF was calculated to be at 2.4. The observed PDF 
is found to be strongly bimodal with peaks at 1.8 and 3.4. 
The strong bimodality is attributed to the very narrow band 
incident wave system and the generation of secondary waves 


Beene  birst Nanmonic. 


Cae Orih PROBABTLILY DENSITY BUNCTIONS (T vs 4H) 

The joint PDF of periods and heights for the plunging 
and spilling examples are contoured in Figures 5 and 6; the 
distributions for all cases are found in Appendix C. 

It is noted that the longest wave periods are associated 
with ene greatest wave heights suggesting strong correlation 
between periods and heights. The correlation coefficient 
was calculated for all cases and varied between 0.60 and 0.80. 
Bretschneider (1959) examined the two special cases of joint 
Rayleigh distributions with correlations of 0 and +l. Per- 
Fect correlation between two variables would result in 
Eomrelarion cOctELcient of 1.0 and a straight line in the 
jem. LDF . 

The joint Rayleigh PDF for the case of a correlation co- 
Peebebeonneon cero 1s plotted im Figure 3. Zero correlation, 
or statistical independence, results in the joint PDF being 
Spread out as compared with the straight line for perfect 
correlation. The amount of spread is a measure of the reduc- 
Eton, im correlation. 

Figure 5 is an example of p(H,T) for plunging breakers 


and shows two regions of high density. The bimodality of 
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p(H) and p(T) are reflected in the p(H,T). The maximum peak 
period, which corresponds to the peak in the wave spectrum, 
is highly correlated with the primary peak of wave heights. 
The second high density region corresponds to the first 
harmonic of the primary wave period. The joint PDF suggests 
strong wave-wave interactions in a plunging wave system re- 
sulting in energy transfer to harmonics. 

Two distinct contour groups are no longer apparent in 
Figure 6 for the spilling case, but rather one long contour 
is found. Since the height PDF and the period PDF are only 
weakly bimodal, it follows that a single contour group 
should be found. The area of maximum joint density is found 
to be associated with the first harmonic of the offshore 
wave period, and the mode of the observed breaker heights. 

Greater density was found at short periods and heights 
for all records than predicted by the Rayleigh PDF. A high 
frequency wave riding on a long period wave causes an addi- 
tional zero crossing as the longer wave passes through zero. 
Most previous comparisons of heights or periods of waves have 
been done using pressure sensors which greatly filter the 
high frequencies due to hydrodynamic attenuation. Capaci- 
tance wave staffs and electromagnetic flowmeters have good 
response times and thus measure high frequencies. Therefore, 
even though the data were filtered at 1 Hz, all records showed 


peaks at the shortest periods and heights. 
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Dee veLeGily PROBABIMIGY DENSITY FUNCTIONS 

The velocity PDFs were calculated as the maximum onshore 
velocities within the period segments determined by the zero- 
up-crosses of the surface elevation record as shown in Fig- 
ures 8 and 9, Using this technique can result in both posi- 
tive and negative maximas for onshore particle velocities 
(i.e., d@v/dt < 0). 

The velocity PDFs for spilling breakers and plunging 
breakers are shown in Figures 8 and 9. The theoretical 
Rayleigh PDF is superimposed on the same graphs for compari- 
son. It is apparent that neither the spilling or plunging 
case is adequately described by the Rayleigh PDF. The veloc- 
ity PDF for spilling breakers appears to be symmetrical, with 
the mode value falling at the bin with center velocity of 
1.5. There also appears to be a small peak value at 0.3. 
When compared with the plunging case, it is found that the 
peaks occur at the same normalized velocities. The plunging 
wave PDF is bimodal and has relatively few occurrences of 


negative onshore velocities. 


Eoe JOINT, PROBABILITY DENSITY FUNCTIONS (V¥ ws 4H) 

The joint probability distributions of velocities and 
heights are shown in Figures @ and 9. The results show high 
correlation between velocities and heights. The high posi- 
tive correlation suggests that the larger maximum onshore 
velocities are correlated with the larger wave heights. 

Figure 8 shows the plunging breaker case. Two peak areas 


ase Golmd £O occur Im the joint PDF which reflects the 
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Dumodailicy Of the indivadual PPPs. The joint PDF for the 
plunging waves suggests a highly correlated case (0.68) in 
which the largest onshore velocities occur due to waves of 
greatest height. 

The spilling case as shown in Figure 9 also shows high 
correlation (0.70). For this example, there is only one 
major peak which corresponds to the individual PDFs which 
are essentially wnimodal. The peak area of the joint PDF 
has been shifted toward higher waves and greater velocities 


than the plunging case. 


foe JOUNT FROBABILITY DENSITY FUNCTIONS (7 vs V) 

Figures 10 and 11 show joint PDFs of periods and veloc- 
ities in which the periods were determined from the zero-up- 
crosses of the wave gauge record and the coincident maximum 
onshore velocities are determined from the flowmeter record. 
The results are similar to the joint PDF of period and 
heights. Spilling waves, as contoured in Figure 11, show one 
Main contour group. The maximum spilling probability is 
found at the intersection of period 0.84 and velocity 1.5. 
The plunging waves as found in Figure 10 reveal two contour 
peaks. One large contour group is associated with the point 
representing the offshore period 1.56 and a velocity of 1.5. 
The second large contour is a function of the first harmonic 
wave period of 0.84 and a velocity of 0.3. 

The negative values of the maximum velocity indicate 
that, during a wave period defined by the zero-up-crosses 


of the surface elevation, the maximum velocity can be negative. 
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The joint PDF of wave height and velocity shows that the 
negative velocities are associated with the low wave heights. 
The joint PDF of velocity and period shows that small nega- 
tive velocities are correlated with the short periods. 
Therefore, it is concluded the negative velocity maxima 

are associated with short periods and low wave heights. The 
secondary waves occur after the passage of the primary wave 
crest and have a peak near the trough of the primary wave; 
the velocities associated with the primary wave would be 
negative maximum in this region so that the measured veloc- 
ity, which is a sum of the primary and secondary waves, can 


be small or negative. 
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V. CONCLUSIONS 


Waves in the surf zone are a highly non-linear process 
which is evidenced by the appearance of secondary waves. 
The secondary waves appeared in the PDFs at the first har- 
monic of the primary period. 

The joint probability density functions of the breaking 
wave properties show high correlations. The correlation 
factor ranged between 0.60-0.80 for periods and heights 
and suggests that the greater the wave period, the higher 
the breaker height will be. High correlation between the 
periods and velocities and the heights and velocities sug- 
gests that the maximum onshore velocities are a function 
of wave periods and breaker heights. 

Breaking waves may be viewed as a gradation between 
two extremes, spilling and plunging breakers. The differ- 
ence between the spilling and plunging breakers becomes 
apparent in the joint PDFs. The spilling breakers yield a 
joint PDF which has only one high density contour region, 
whereas the plunging breakers show two high density contour 
groups. The high density contour groups appear to be most 
associated with the period PDF. 

The joint PDF (p(T,H)) of spilling breakers has a maxi- 
Mum density at a wave period equal to one half the peak 


period of the wave energy density spectrum, i.e. the first 
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harmonic. The corresponding velocities (p(v,T)) or heights 
(p(H,T)) were likewise found to be unimodal with a peak den- 
sity which is driven by the first harmonic wave period. 

Plunging breakers, at the opposite end of the gradation, 
yield a bimodal joint PDF. The two major. contour groups 
appear to be a result of the bimodal density distribution 
of the periods. The two peak densities of the period PDF 
are representative of the peak of the wave density spectrum 
and a first harmonic. Due to the high correlation of periods 
and heights, the corresponding PDF of wave heights appears 
to be bimodal. The joint PDF of period and velocity and the 
joint PDF of velocity and height suggests high correlation 
between these variables. The maximum onshore velocities 
appear to be correlated with both wave periods and wave 
heights. 

A truncated Rayleigh PDF was found to loosely approxi- 
mate the empirical period and the height PDFs of spilling 
breakers. This truncation seems to be the result of long 
period waves with their associated larger breaker heights 
breaking prior to the arrival at the wave gauge. The mode 
of the spilling wave period distributions corresponds to the 
first harmonic of the peak of the energy density spectrum. 
The smoothed appearance of the spilling wave periods PDF, 
relative to the plunging case, appears due to the smearing 
which occurs with the appearance of harmonics greater than 
the first. The velocities PDFs of spilling waves, where the 


velocity maximum is determined as in Section D, yields a 
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symmetrical distribution with some negative onshore veloc- 
ities present. Since the Rayleigh distribution is undefined 
for negative values, the resulting fit of observed data with 
the theoretical Rayleigh distribution is poor. 

The PDFs of theplunging breakers are poorly described 
by either the Rayleigh PDF or other existing theoretical 
PDFs such as the Weibull PDF. Theoretical PDFs are unimodal, 
whereas the empirical PDFs under plunging conditions are 
strongly bimodal. The high correlation found in the joint 
PDFs for wave heights and maximum velocities also results 
in a bimodal distribution both in the velocities and the 
heights PDFs. 

The mean wave period as obtained by the zero-up-cross 
method in the surf zone is a poor indicator of the mean 
offshore wave period. Spectral analysis of the breaking 
waves shows an energy peak corresponding to the primary off- 
shore wave period, but the period PDFs of breaking waves 
have a strong peak at the first harmonic of the wave spectra 
peak. Plunging waves had high density peaks in the period 
PDF corresponding both to the peak of the wave spectrum and 
to its first harmonic. The result is that the mean of the 
breaking waves falls between the primary and the first har- 
monic wave period. The mean period of breaking waves there- 
fore is a poor estimation of the mean period of the incident 


Wavy eS. 
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APPENDIX 


STATISTICAL PARAMETERS CALCULATED 


1. Due to the variation of sampling frequencies between the 
different experiments, it follows that some of the other 
statistical parameters will also vary. A summary of these 
parameters follows where fn is the Nyquist frequency, Af 
is the resolution, At is the interval between samples and 


T ais the length of sampled record in minutes. 





DATE fy (iz) Af (Hz) At(sample/sec) T(min) 
4 MARCH 1975 2.5 £6.1x107* 0.200 Nee 
5 MARCH 1975 2.5  6.1x107* 0.200 Ape 
PeMRCKTo 75) 9255  6.1x107* 0.200 D7 
8 MARCH 1975 2.0 4.8xl07" 0.256 35.0 

29 MAY 1975 2.0 4.9x107% 0.256 35.0 

16 MARCH 1977. 2.0 4.9x107% 0.250 ae 

17 MARCH 1977. 2.0 4.9xl07* 0.250 aol 

19 MARCH 1977. 2.0 4.9xl0~¢ 0.250 Vs ik 
8 MARCH 1978 2.0 4.9x107* 0.250 gael 
9 MARCH 1978 2.0 4.9x107° 0.250 ek 


2. Tables of calculated statistical parameters follow where 


Nn = sea surface elevation (m) 

v = instantaneous particle velocities (m/sec) 

VM = maximum onshore particle velocities (m/sec) 

LR = length of record for calculated statistics (min) 
—€ = spectral width parameter 

H = wave height (m) 

T = wave period (sec) 
ig 


= correlation coefficient 


at 
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APPENDIX B 


MATHEMATICAL DERIVATION OF NORMALIZED RAYLEIGH DISTRIBUTIONS 


Wave heights, periods and velocities have been compared 


with the Rayleigh distribution 


2 
A A 
Sar <P ala) 
46 26 
where A is the non-normalized value. Since all parameters 


are normalized in this paper, the following conversion was 
made: 


let A' = A/o , 














then dA . 
p(A') = p(A) (Gi) A=A'o 
p(A') = p(a) S42 = prayo 
for wave heights 
p(H) = +, exp (-4(s=-)*) 
4o n 
n 
H. 4H Heo 
p(—)= pe exp (-k(g=)°) 
ia) ry ia) 
for maximum velocities 
V 2 
p(V) = =) exp(-(s+)*) 
Lo Vv 
V 
V aw 
ree yee ee) ) 
 oehe - 2o., 
for periods 
4g 
p(T) = ie exp (-%(sx)*) 
O 


54 





where 


then 


li 


rye 


DS 








APPENDIX C 


PROpabiE ity DENSLTY BUNCTIONS 


The following tables consist of both the normalized 
joint probability density functions, and the individual 
probability density functions as leipanlee es The PDF tables 
coincide with Figures 5 thru ll, respectively. 

The top table is the joint PDF. The top line and the 
left column represent the center period, height, or 
velocity point of each of the 15 bins. The numbers in 
the other columns and lines represent the frequency of 
occurrence in the respective bins. 

The bottom table represents the individual PDFs as 
labeled. The top line of the individual PDFs is the 
center point of the individual bins. The second line 


represents the frequency of occurrence in the bin. 
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